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It is conventional to choose a typical momentum transfer of the process as the renormalization 
scale and take an arbitrary range to estimate the uncertainty in the QCD prediction. However, 
predictions using this procedure depend on the renormalization scheme, leave a non-convergent 
renormalon perturbative series, and moreover, one obtains incorrect results when applied to QED 
processes. In contrast, if one fixes the renormalization scale using the Principle of Maximum 
Conformality (PMC), all non-conformal {,3;}-terms in the perturbative expansion series are summed 
into the running coupling, and one obtains a unique, scale-fixed, scheme-independent prediction 
at any finite order. The PMC scale /i R and the resulting finite-order PMC prediction are both 
to high accuracy independent of the choice of initial renormalization scale fJ, 1 ^ 11 , consistent with 
renormalization group invariance. As an application, we apply the PMC procedure to obtain NNLO 
predictions for the it-pair production at the Tevatron and LHC colliders. The PMC prediction for 
the total cross-section a t i agrees well with the present Tevatron and LHC data. We also verify that 
the initial scale-independence of the PMC prediction is satisfied to high accuracy at the NNLO 
level: the total cross-section remains almost unchanged even when taking very disparate initial 
scales fi 1 ^ equal to m t , 20 m t , ^/s. Moreover, after PMC scale setting, we obtain A l p B ~ 12.5%, 
A p / B ~ 8.28% and Af B (M tl > 450 GeV) ~ 35.0%. These predictions have a 1 cr-deviation from 
the present CDF and DO measurements; the large discrepancy of the top quark forward-backward 
asymmetry between the Standard Model estimate and the data are thus greatly reduced. 
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Physical predictions in Quantum Chromodynamics 
(QCD) are in principle invariant under any choice of 
renormalization scale and renormalization scheme. It is 
common practice to simply guess a renormalization scale 
H>R — Q, Q being a typical momentum transfer of the 
process, and then vary it over the range [Q/2, 2 Q]. How- 
ever, this arbitrary procedure leads to scheme-dependent 
predictions at any finite order in perturbation theory. In 
fact, a principal ambiguity in perturbative QCD calcula- 
tions lies in the choice of fiR. It has been considered as 
a main systematic error in QCD perturbative analyses. 

The Brodsky-Lepage-Mackenzie method (BLM) 1| 
and the Principle of Maximum Conformality (PMC) [2j,|3| 
provide a solution to this problem. The PMC provides 
the principle underlying BLM scale setting; the BLM is 
equivalent to PMC through the PMC - BLM correspon- 
dence principle [3] , so we shall treat them on equal foot- 
ing. When one applies the PMC, all non-conformal {ft}- 
terms in the perturbative expansion are summed into the 
running coupling so that the remaining terms in the per- 
turbative series are identical to that of a conformal the- 
ory, i.e., the corresponding theory with {Pi} = {0}. 

The PMC coefficients and PMC scales may be dif- 
ferent under different renormalization schemes, however 
their combined result will be the same, since the scheme- 
dependent PMC scales for different schemes are related 
by commensurate scale relations J4j. Thus, QCD predic- 
tions using PMC are independent of the renormalization 
scheme. After PMC scale setting, the divergent "renor- 



malon" series with n!-growth disappear, so that a more 
convergent perturbative series is obtained. 

The PMC method satisfies all self-consistency con- 
ditions, including the existence and uniqueness of the 
scale, reflexivity, symmetry and transitivity [5]. In the 
Abelian limit Nc — > at fixed a = Cf&s with Cf — 
[Nl - 1)/2N C , it agrees with the Cell Mann-Low pro- 
cedure for setting the scale in QED [y, |7|. Thus as in 
QED, the renormalization scale can be unambiguously 
set at each finite order by the PMC. The PMC scales 
and coefficients can be set order-by-order. A systematic, 
scheme-independent procedure for setting PMC scales up 
to next-to-next-to-leading-order (NNLO) has been pre- 
sented in Ref. J3J. 

Formally, one needs to choose an initial renormaliza- 
tion scale /U™* for PMC. However, the final result when 
summing all {/3i}-terms to all orders will be independent 
of /xg H ; i.e. for any observable O, dO (^ MC ) / d ^ = 0, 
where /^ MC stands for the PMC scale. This is the in- 
variance principle used to derive renormalization group 
results such as the Callan-Symanzik equations [8]. The 
PMC scales in higher orders take the form of a pertur- 
bative series in a s so as to properly absorb all {/?i}— de- 
pendent terms associated with renormalization into the 
a s -running coupling [3J, |j]. At fixed order, there is some 
residual initial-scale dependence because of the unknown- 
higher-order {/^}-terms. However, such residual renor- 
malization scale-uncertainty will be greatly suppressed 
since those higher order {/3i}-terms will be absorbed into 



the PMC scales' higher-order a s -terms. 

As an important application of the PMC, we shall 
predict the it-pair hadroproduction cross-section a t i up 
to NNLO. It has been measured at the Tevatron and 
LHC with high precision |9T4l2J[. Theoretically, o t i has 
been calculated up to NLO within the MS'-scheme 13 1. 
Large logarithmic corrections associated with the soft 
gluon emission have been investigated and resumme d up 
to next-to-next-to- leading- logarithmic corrections 14J . 



Even though complete NNLO fixed-order results are not 
available, parts of the fixed-order NNLO results have 
been derived through resummation 15]. These results 
provide the foundation for estimating the NNLO results. 
The hadronic cross-section for top quark pair produc- 
tion can be written as: 
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where the parton luminosity 



lij = S ~^ /Hl fc 1 '^) /?'/#= ( X2 'W) 



with X\ = s/S and x 2 — s/s. Here S denotes the 
hadronic center-of-mass energy squared and s = X1X2S 
is the subprocess center-of-mass energy squared. The 
parameters pu and Pf denote the renormalization and 
factorization scales, and the functions f i / Hl2 (x a ,pf) 
(a = 1,2) are the parton distribution functions (PDFs) 
describing the probability to find a parton of type i with 
a momentum fraction between x a and x a + dx a in the 
hadron H\.i- The top quark mass rat is the mass renor- 
malized in the on-shell scheme. 

The partonic subprocess cross-sections &ij can be de- 
composed in terms of the dimensionless scaling-functions 
/™, where (ij) = {{qq),(gg),{gq),{gq)} stands for the 
four production channels and m = 0, 1, 2 stands for LO, 
NLO and NNLO functions respectively. The analytical 
expressions for f t J ' (ft, Q) which contain the explicit fac- 
torization and renormalization scale dependence can be 
directly read from the HATHOR program [16j. Up to 
NNLO, &ij takes the following form 



^•=^i;/^.o)rt) 



(2) 



where p = 4ra^/s and a s (Q) = a s (Q)/ir. There is un- 
certainty in setting the factorization scale pf which ap- 
pears even in conformal theory, and its determination is a 
completely separate issue from the renormalization scale 
setting. To keep our attention on the renormalization 
scale, we implicitly set pf = ra t . For the initial value of 
M-R = Mfl 1 *) we take /i™* = Q, where Q stands for the 
typical momentum transfer of the process. For example, 
Q can be taken as m t , 2ra t , y/s, etc. As the default 
choice, we take Q = mt. 

According to the PMC, we need to identify the 

(1 2) 

n\ ' -dependent terms associated with renormalization. 
Coulomb-type corrections are enhanced by factors of n 
and the PMC scales can b e relatively soft for heavy quark 
velocity v = yl — 4rnf/s — > 0. Thus the terms which are 
proportional to (n/v) or {ir/v) 2 have a separate PMC 
scale and will thus be treated separately [17(. More ex- 
plicitly, the NLO and NNLO scaling-functions can be 
written as 

fl j (p,Q) = [A lij + B lij n f } + D Uj {^) , (3) 

f!j(p, Q) = [A-Xj + B 2lJ n f + C 2lJ n)] + 

[D 2ij +E Mj n f }(^)+F 2ij (~y. (4) 

The PMC scales can be set order-by-order and the final 
result is 
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where k = ^a s {Q* 2 ) + ^f- (f ) a s {Q* 2 ). Here Q\ and 
Ql* are the LO and NLO^PMC scales for the non- 
Coulomb part, and Q 2 is the LO PMC scale for the 
Coulomb part. The PMC coefficients and PMC scales, 
together with their detailed derivations, can be found in 
Ref. [dj. 

When we do numerical calculations, the input param- 
eters are chosen with the following values: for the top 
quark mass, we adopt the PDG value [19(; i.e. m t — 
172.9 ± 1.1 GeV. For the PDFs, we adopt the CTEQ 
CT10 set with a s {m z ) = 0.118 [H. Our results for 
the ti production cross-sections are presented in Table U 
where the total cross-sections which are derived by using 
the PMC scale setting and the conventional scale setting 
are presented. 



1 m=0 



It is found that after PMC scale setting, the result- ing total cross-sections for five disparate initial scales are 





PMC scale setting 


Conventional scale setting 




Q - m t /4 


Q = m t 


Q = 10m t 


Q = 20m t 


Q=V^ 


[ir = mt/2 


/in = m t 


fi R = 2 m t 


Tevatron (1.96 TeV) 


7.620(5) 


7.626(3) 


7.625(5) 


7.624(6) 


7.628(5) 


7.742(5) 


7.489(3) 


7.199(5) 


LHC (7 TeV) 


171.6(1) 


171.8(1) 


171.7(1) 


171.7(1) 


171.7(1) 


168.8(1) 


164.6(1) 


157.5(1) 


LHC (14 TeV) 


941.8(8) 


941.3(5) 


942.0(8) 


941.4(8) 


942.2(8) 


923.8(7) 


907.4(4) 


870.9(6) 



TABLE I. Dependence of the tt production cross-sections (in unit: pb) at the Tevatron and LHC on the initial renormalization 
scale fj^R 1 — Q. Here m t — 172.9 GeV. The number in parenthesis shows the Monte Carlo uncertainty in the last digit. 
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FIG. 1. Total cross-section a t i for the top quark pair produc- 
tion versus top quark mass. 



equal to each other within part per mill accuracy q For 
comparison, we also present the results with conventional 
scale setting in Table [U For \ir e [m t /2,2m t ], we ob- 
tain the usual renormalization scale-uncertainty I _ 4 y ) . 
This shows that the renormalization scale uncertainty is 
greatly suppressed and essentially eliminated using PMC 
even at the NNLO level. This is consistent with renor- 
malization group invariance: there should be no depen- 
dence of the prediction for a physical observable on the 
choice of the initial renormalization scale. 

The PMC predictions for total cross-section <j t t are 
sensitive to the top quark mass. We present a t i as a 
function of mt in Fig. ((J}. After PMC scale setting, the 
value of a t t becomes v ery close to the central values of the 
experimental data [9h12||. By varying mt = 172.9 ± 1.1 
GeV [ljj, we predict 
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FIG. 2. Comparison of the PMC prediction with the CDF data [21( for the ii-pair forward-backward asymmetry for the whole 
phase-space. The left diagram is for A f p B in the ii-rest frame, the middle diagram is for A V p B in the laboratory frame, and 
the right diagram is for A t p B (M t t > 450 GeV). The Hollik and Pagani's results (HP) 24] using conventional scale setting are 
presented for a comparison. The result for DO data [23] shows a similar behavior. 



1 There is some small residual initial-scale dependence in the PMC 



scales because of unknown-higher-order {/3j}-terms. 
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We have recently shown that the large discrepancy be- 
tween the Standard Model estimates using conventional 
scale setting and the CDF and DO data 21. 22] for the 
it-pair forward-backward asymmetry is mainly caused 
by improper setting of renormalization scale [23]. Af- 
ter PMC scale setting, it is found that the NLO PMC 
scale has a dip behavior for the dominant asymmet- 
ric (gg)-channel; the importance of this channel to the 
asymmetry is thus increased. Then, after PMC scale 
setting, the it-pair forward-backward asymmetries Ap B 
and A P p B at the Tevatron are increased by 42% in com- 
parison with the previous estimates obtained by using 
conventional scale setting. We obtain Ap B ~ 12.5%, 
A% ~ 8.28% and A% B {M tl > 450 GeV) ~ 35.0% Q. 
These predictions have a la-deviation from the present 
CDF and DO measurements; the large discrepancy of 
the top quark forward-backward asymmetry between the 
Standard Model estimate and the data are thus greatly 
reduced. This large improvement is explicitly shown in 
Fig-©i where Hollik and Pagani's results, which are de- 
rived under conventional scale setting 24| , are presented 
for comparison. 

Summary: By using PMC scale setting, one obtains 
a unique, scale-fixed, scheme-independent prediction at 
any finite order in a systematic way. Since the renormal- 
ization scale and scheme ambiguities are removed, this 
procedure improves the precision of tests of the Stan- 
dard Model and enhances the sensitivity to new phe- 
nomena. The PMC can be applied to a wide-variety of 
perturbatively-calculable collider and other processes. 

We have applied PMC to study the tt hadroproduc- 
tion cross-section a t i up to NNLO. The resulting LO- 
and NLO- terms are conformally invariant and scheme- 
independent, and the non-conformal contributions in the 
NNLO-terms are greatly suppressed. The PMC predic- 
tion for a t t agrees well with the present Tevatron and 
LHC data. We also verify that the initial renormalization 
scale-independence of the PMC prediction is satisfied to 
high accuracy at the NNLO: the total cross-section re- 
mains almost unchanged even when taking very disparate 
initial scales /z™* equal to m t , 10 m t , 20 m t , y/s. The 
optimized PMC scales substantially eliminates the large 
discrepancy between the Standard Model estimation and 
the Tevatron data for the it-pair forward-backward asym- 
metry. 
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